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Abstract. Worst-case flutter margins may be computed for 
a linear model with respect to a set of uncertainty operators 
using the structured singular value. This paper considers an 
on-line implementation to compute these robust margins in a 
flight test program. Uncertainty descriptions are updated at 
test points to account for unmodeled time- varying dynamics 
of the airplane by ensuring the robust model is not invali- 
dated by measured flight data. Robust margins computed 
with respect to this uncertainty remain conservative to the 
changing dynamics throughout the flight. A simulation clear- 
ly demonstrates this method can improve the efficiency of 
flight testing by accurately predicting the flutter margin to 
improve safety while reducing the necessary flight time. 
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1. Introduction 

Flight flutter testing incurs dramatic time and 
costs associated with safely expanding the flight 
envelope to ensure no aeroelastic instabilities are 
encountered. NASA Dryden Flight Research Cen- 
ter is investigating methods to compute confident 
flutter margins to improve flight test efficiency [11]. 

Traditional pre-flight analytical methods of pre- 
dicting flutter margins, such as the p-k method, 
utilize a single model of the aircraft dynamics. 
These methods are not able to account for vari- 
ations between the model and the aircraft so the 
predicted flutter margins may be overly optimistic. 

Traditional flight test methods of determining flut- 
ter margins often utilize tracking estimates of 
modal damping obtained from flight data [6]. These 
methods are both inefficient and dangerous. Mul- 
tiple data sets must be taken at each test point to 
ensure the critical modal dynamics have been suffi- 
ciently excited. Also, the envelope must be expand- 
ed in small increments since damping only indi- 
cates stability at the current test point and cannot 
be extrapolated to consider other flight conditions. 

Several proposed on-line methods use modal filter- 
ing, parameter identification and envelope analysis 
to increase confidence in the flutter boundary [3]. 
These methods, like damping tracking, compute 
stability indicators rather than predictors. 


A novel method for computing flutter margins is 
developed based on robust stability theory [8]. 
Uncertainty operators are associated with a linear 
model to account for modeling errors. The struc- 
tured singular value, /x, computes a stability mar- 
gin robust to these uncertainties. In this respect, a 
worst-case flutter margin is computed. 

This paper considers an on-line implementation 
using /x, a stability predictor, to replace damp- 
ing. This method presents advantages over tradi- 
tional flutter analysis methods since it combines 
a well developed analytical model with flight data 
describing the true aircraft, as shown in Figure 1 . 
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Figure 1. Methods for Computing Flutter Margins 

A simulated flight test demonstrates several signif- 
icant advantages of the fi method as compared to 
traditional methods. The /x method is able to gen- 
erate realistic flutter margins using flight data from 
a stable flight condition; whereas little information 
is obtained from tracking damping until the onset 
of flutter. Also, the robust flutter margins are con- 
sistently conservative due to model validation of 
the uncertainty set using flight data; whereas nom- 
inal margins computed with p-k do not account for 
changes in aircraft dynamics. 

The concept of a “flutterometer* is introduced as 
a flight test tool to indicate the proximity of a flut- 
ter condition. The simulation demonstrates a fea- 
sible concept of such a tool based on robust flutter 
margin predictions which remain conservative to 
time-varying aircraft dynamics. 
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2. Robustness Measure /z 


3, Robust Flutter Model 


Any aeroelastic model is an approximate repre- 
sentation of the aircraft dynamics. Uncertainty 
operators are associated with the model in the 
framework of the structured singular value, /z, to 
account for inaccuracies such as errors and unmod- 
eled dynamics. 


Consider the aeroelastic equation of motion [5], 

Mr) + Cri + KTi + qQ{s)r} = 0 (1) 

where the unsteady aerodynamic forces are 

Q(s) = D Q + C q (si -Aq)- 1 Bq (2) 


The uncertainty operator, A, is allowed to lie with- 
in a norm bounded set. Weighting matrices are usu- 
ally included in the /z framework to normalize the 
uncertainty norm bound to unity. 

A = {A : ||A||oo < 1} 

Define /z with respect to the set A. 

^ = min {ct(A) : A € A, det(/ - PA) = 0} 

with /z(P) = 0 if there exists no A 6 A such that 
det(/-PA) =0. 

Define elements of P as transfer function from dis- 
turbances to errors P 22 , transfer function relating 
uncertainty feedback Pn, and similar P 12 and P 21 . 

/z is an exact measure of robustness for systems 
with structured uncertainty. The inverse of /z can 
be interpreted as a measure of the smallest destabi- 
lizing perturbation. The system is guaranteed to be 
robustly stable for all uncertainty operators bound- 
ed by the smallest destabilizing value [4]. 

Theorem 2.1. A stable system P is robustly stable 
for all A E A if and only if /z(Pn) < 1. 

Model validation algorithms can verify the uncer- 
tainty description is sufficient to describe varia- 
tions between the model and measured data. A 
/z- analysis method of model validation has been 
developed that determines if there exists some 
perturbation to the system dynamics, within the 
allowed set of uncertainty, for which the robust sys- 
tem could produce the measurement data for the 
indicated input data [7], 

Theorem 2.2. Given measurements y generated by 
input tz, then a system P with uncertainty A is not 
invalidated if /z (P\\ — P 12 U (P 22 U “ y) 1 P 21 ) > 1 


Unfortunately, /z is difficult to compute. Upper 
and lower bounds have been derived which may 
be solved efficiently using convex optimization and 
a power iteration [1], 


This system is parameterized around q such that 
the /z value describes the smallest variation in 
dynamic pressure for the onset of flutter. This is 
accomplished by considering an additive perturba- 
tion, $g, on the nominal dynamic pressure q Q . 

9 = 9o + % (3) 


The nominal plant P n0 m , based on Equations 1 
and 2, incorporates the perturbation ^ as a lin- 
ear fractional transformation using signals z and w 
where w — &qZ is the perturbation feedback rela- 
tionship. Stability with respect to a range of 8q val- 
ues is analyzed with the robust stability condition 

of Theorem 2.1 with Pnjwl, such that z — PnJlw, 
is given below. Define M — —M~ l . 
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The full nominal plant, P n0 m, utilizes Pn^ with 
additional elements to describe the generation of 
measurements and effects of force excitation. 


A robust model P ro f, associates uncertainty oper- 
ators with Pnom to account for modeling errors 
between the theoretical system and the physical 
aircraft. These operators may enter as paramet- 
ric uncertainties to reflect variations in specific 
elements such as mass and damping or dynamic 
uncertainties to allow unstructured variations in 
magnitude and phase of signals. 

Aeroelastic flight data, if available, can be used 
to determine a realistic uncertainty description for 
P ro * using Theorem 2.2. Utilizing flight data may 
increase confidence in the robust model; however, 
it may also introduce additional uncertainty due to 
flight test issues, aircraft performance, aeroelastic 
excitation and frequency domain mapping [11, 9]. 

Robust flutter margins are computed for P ro b using 
the robust stability condition of Theorem 2.1. A 
robust flutter margin is the smallest destabilizing 
perturbation, &q, to dynamic pressure for P ro & with 
the given amount of uncertainty. This is the worst- 
case margin for the range of considered dynamics. 
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4. Robust Flutter Parameter 

The flutter computation method described in this 
paper uses /x as the worst-case flutter parameter. 
There are several advantages to using fi as the flut- 
ter parameter, fi is a much more informative flutter 
margin as compared to traditional parameters such 
as pole location and modal damping. 

The conservatism introduced by considering the 
worst-case uncertainty perturbation can be inter- 
preted as a measure of sensitivity. Robust /x values 
which are significantly different than the nominal 
flutter margins indicate the plant is highly sensitive 
to modeling errors and changes in flight condition. 
A small perturbation to the system can drastical- 
ly alter the flutter stability properties. Conversely, 
similarity between the robust and nominal flutter 
margins indicates the aircraft is not highly sensi- 
tive to small perturbations. 

Robustness analysis determines not only the norm 
of the smallest destabilizing perturbation but also 
the direction. This information relates exact per- 
turbations for which the system is particularly sen- 
sitive. fi can thus indicate the worst-case flutter 
mechanism which may naturally extend to indicate 
active and passive control strategies for flutter sup- 
pression. 

Additionally, damping is only truly informative 
at the point of instability since stable damping 
and damping trends at a given flight condition do 
not guarantee what increases in dynamic pressure 
may be safely considered, /i computes the small- 
est destabilizing perturbation which indicates the 
nearest flight conditions that will cause a flutter 
instability. In this respect, /x is a stability predic- 
tor while damping is merely a stability indicator. 

These characteristics of fi make the worst-case flut- 
ter algorithm especially valuable for flight test pro- 
grams. Aeroelastic flight data can be measured 
at a stable flight condition and used to evalu- 
ate uncertainty operators. The /x method, unlike 
damping estimation, does not require the aircraft 
to approach instability for accurate prediction, /x 
can be computed to update the stability margins 
with respect to the new uncertainty levels. The 
worst-case stability margin then indicates what 
flight conditions may be safely considered for safe 
and efficient expansion of the flight envelope. 


5* On-Line Implementation 

Traditional flight flutter testing for envelope clear- 
ance proceeds with several repeated operations [6]. 
The aircraft is trimmed at a stable flight condition 
on the edge of the envelope where measured aeroe- 
lastic responses are analyzed to estimate damping 
and modal characteristics. If the estimates do not 
indicate the imminent onset of flutter, the enve- 
lope is slightly expanded and the procedure repeats 
until flutter or a performance limit is reached. 

Worst-case flutter margins can easily be computed 
on-line in such a flight test program. The imple- 
mentation can compute these margins in parallel 
with the damping estimates. Both methods can uti- 
lize the same data to update the stability measure 
at each test point. The robust model, which uses 
the same nominal plant throughout the flight, sim- 
ply updates the uncertainty set such that the model 
is not invalidated by the flight data and then com- 
putes the worst-case margin with respect to that 
uncertainty. 

Several robust flutter margins may be computed 
at a test point by computing margins which are 
worst-case to different uncertainty sets. Two obvi- 
ous approaches are denoted here as local or global. 

A local approach computes an uncertainty set at 
each test point such that the model is not invalidat- 
ed by the current flight data. No knowledge of past 
uncertainties or data is used. This approach may 
seem logical considering a linear transonic mod- 
el, which is often suspect, associates an uncertain- 
ty description which may be overly conservative 
for a subsonic model. This method may be highly 
accurate but could be overly optimistic if critical 
dynamics are not sufficiently excited. 

A global approach computes flutter margins with 
respect to the largest uncertainty required for any 
point in the flight. It may be reasonable that if a 
model has a certain level of error at a given flight 
condition, then that same level might exist at all 
flight conditions. This method may be more con- 
servative but the safety of the margin is increased 
since it is less susceptible to poor excitation at a 
certain test point. 

Many flight analysis facilities are increasing com- 
putational resources to decrease required flight 
time. Such environments make on-line analysis of 
robust flutter margins extremely attractive. Sepa- 
rate computers can be used so worst-case flutter 
margins do not replace damping estimates; rather, 
they augment traditional flutter measures. 
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6. Simulation : Aircraft Model 

A supersonic flight flutter test is simulated to 
demonstrate the benefits of utilizing the robust 
flutter margin algorithm in an on-line manner. 

The F/A-18 Systems Research Aircraft (SRA) pro- 
vides the aircraft model for the simulated flight 
test. This aircraft is a two-seat configuration fight- 
er with production engines being flown at NASA 
Dryden Flight Research Center [2]. 

There are three linear state-space models used in 
the simulation. 

Pnom • nominal theoretical model 
P TO b : robust theoretical model 
Ptrue : true aircraft dynamics 

The nominal theoretical model, P na m , is a struc- 
tural model coupled with unsteady aerodynam- 
ic forces. Pnom is parameterized around dynamic 
pressure using an operator ^ as in Equation 4. 
This model was used for flutter analysis of the F/A- 
18 SRA using the n method and traditional algo- 
rithms such as the p-k method [8, 10, 12], 

Pnom is a 84 state plant for this simulation to model 
the 14 symmetric structural modes of the aircraft 
along with the 56 states used in the state-space 
representation of the unsteady aerodynamic forces. 

The robust theoretical model, P ro b, associates an 
uncertainty description with the nominal plant 
to describe errors and unmodeled dynamics. Two 
uncertainty operators, A a and A are included in 
the robust model to account for variations between 
the nominal plant and the true aircraft dynamics. 

The uncertainty operator A^ affects the state 
matrix of the nominal plant to model variations 
in both natural frequency and damping for each 
mode. A^ is a struct ured diagonal matrix with real 
scalar parameters as elements. Separate scalars are 
used to affect each modal response and time lag in 
the state matrix. A scalar associated with a modal 
response is repeated two times while each time lag 
uncertainty appears once on the diagonal. 

Initial weights for these uncertainties are chosen 
based on previous analysis of P n0 m at superson- 
ic flight conditions [8, 10]. The weighting W u 
describes the percent of variation allowed in natu- 
ral frequency while W$ relates to damping and Wi 
is associated with time lags in the unsteady aero- 
dynamic forces. 

Wu, = .05 W c = .15 W { = .15 (5) 


The uncertainty operator A i n is a complex mul- 
tiplicative uncertainty on the excitation force. A 
weighting function, W{ n , reflects the frequency 
varying levels of multiplicative uncertainty such as 
a large component at high frequency to indicate no 
dynamics above 40 Hz are included in the model. 


Win = 5 


s + 100 
s + 5000 


( 6 ) 


P ro f, also contains sensor noise with a magnitude 
up to 10% of the measurement signal. 


The block diagram for robust flutter margin anal- 
ysis of the F/A-18 SRA model P TO b, including the 
parametric variation in dynamic pressure &q and 
uncertainties A^ and Aj n , is given in Figure 2. 



Figure 2 . Robust Model, P ro &, of the F/A-18 SRA 

The model Ptrue is used as the ‘true’ F/A-18 SRA 
aeroelastic dynamics. Ptrue is based on the nomi- 
nal theoretical plant; however, Pnom has two major 
inaccuracies due to changes in the dynamics of the 
true aircraft. The type and magnitude of the errors 
in Pnom are chosen from within a set of model- 
ing errors determined from post-flight analysis of 
recent flight tests with the F/A-18 SRA [10]. 

— error in structural damping 

— unmodeled time-varying mass 

The structural damping of Ptrue is 10% less than 
that of Pnom - This is introduced by a constant scal- 
ing of 0.9 to the C matrix term of Equation 1. 

The mass of the true aircraft is time- varying 
throughout the simulation while the nominal plant 
is formulated using a constant heavyweight con- 
dition. The dynamics of Ptrue are more realistic 
than Pnom since mass may vary significantly dur- 
ing actual flight testing. 

The Trailing Edge Flap mode for Ptrue at Mach 1,2 
is the critical flutter mode. Flutter conditions vary 
from 27.36 Hz for q = 2360 lb /ft 2 with Ptrue at 95% 
heavyweight mass to 27.38 Hz for q = 2255 lb /ft 2 
at 90% heavyweight mass. 
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7. Simulation : Flight Test 

The simulation is designed to model the procedures 
used in a real flight flutter test. For ease of presen- 
tation, only the simulated flight test to determine 
the flutter boundary at Mach 1.2 for the symmetric 
modes is discussed. Flutter clearance of the entire 
flight envelope including antisymmetric modes is a 
straightforward extension to this simulation. 

The flight test procedure involves iterating over 
several steps. 

1. decrease altitude to raise q by 100 lb /ft 2 

2. measure symmetric response data 

3. determine an uncertainty set such that the 
flight data does not invalidate the model 

4. compute robust flutter margin 

The initial flight condition is q = 100 lb /ft 2 at the 
altitude of 68.8 kft. The aircraft proceeds to the 
first test point at q = 700 lb /ft 2 which corresponds 
to an altitude of 27.7 kft. Further test points occur 
at intervals of 100 lb /ft 2 . 

The aircraft follows the same iterative procedure 
throughout the simulation until the onset of a flut- 
ter instability. Naturally during a real flight test 
the aircraft would not purposely approach a flut- 
ter margin so closely, but encountering the flutter 
condition serves to demonstrate the accuracy of the 
computed flutter margins for the simulation. 

The total flight path is given in Figure 3. 



Figure 3. Flight Path at Constant Mach=1.2 


Time is not a commanded variable in this simula- 
tion so the length of the flight test is determined by 
the computational analysis time. The only known 
function of time is the flight path between test 
points. The aircraft is instructed to increase the 
dynamic pressure by 100 lb /ft 2 and stabilize at 
Mach 1.2 in 30 seconds. 


The horizontal portions of the flight path indicate 
time spent at a test point for which the flight con- 
ditions do not change. The length of time at each 
test point is determined by the computational cost 
of modal validation and robust flutter margin anal- 
ysis. The unequal computational times at different 
test points demonstrated in Figure 3 is due to vari- 
ations in the number of iterations required to gen- 
erate and validate increased uncertainty levels. 

The true aircraft dynamics are time-varying due 
to decreasing mass throughout the simulation. The 
nominal and robust plant models used for flutter 
margin prediction are formulated for the heavy- 
weight flight condition and thus are progressively 
worse representations of the true aircraft dynamics. 
Consequently, the amount of uncertainty associat- 
ed with the robust model must increase to ensure 
the allowed range of dynamic variations includes 
the true dynamics. 

The mass of the aircraft Ptrue is a time-varying 
element of the simulation. The initial mass is cho- 
sen as 95% of the heavyweight condition of the 
SRA with full fuel. The mass variation proceeds 
throughout the simulation by decreasing at a rate 
of 5% of the heavyweight value for 20 minutes 
of flight time. Figure 4 shows the mass of Ptrue 
throughout the simulation. 



Figure 4 . Mass of the Aircraft Ptrue Expressed as Percentage 
of Heavyweight Condition 

The time- varying mass will affect every aeroelastic 
modal response. The decrease in mass will appear 
as increases in natural frequency and damping. The 
modal characteristics for the Trailing Edge Flap 
mode of Ptrue at q— 1000 lb /ft 2 are u = 25.68 Hz 
and £ = .0125 for 95% of heavyweight mass while 
u) = 25.74 Hz and C = 0131 for 90% of heavy- 
weight mass. P TO i, must be robust to these changes 
in modal dynamics and the corresponding change 
in flutter margins of P trlie . 
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The uncertainty levels of P ro & are analyzed at each 
test point to ensure the measured flight data does 
not invalidate the model. If the robust model is 
invalidated, the uncertainty weights are iteratively 
increased and analyzed until the amount of uncer- 
tainty is sufficient to model the observed variations 
between Ptrue and P ro 6- 

This simulation increases the amount of modal 
uncertainty to ensure the flight data could be 
generated by the robust model while keeping the 
amount of input multiplicative uncertainty con- 
stant. The largest increases in uncertainty are for 
the damping levels with smaller increases associat- 
ed with natural frequency uncertainty. The amount 
of uncertainty in damping is expressed as a percent 
of the nominal damping value in Figure 5. 



Figure 5. Level of Modal Damping Uncertainty 


The initial uncertainty levels in the simulation are 
pre-defined based on previous analysis of the F/ A- 
18 SRA aircraft and models. Figure 5 demonstrates 
this initial uncertainty level of approximately 5% 
is sufficient to validate the robust model according 
to flight data recorded between 3 and 8 minutes of 
the simulation. 

Increases in the damping uncertainty are generated 
at several points between the 8 t/l and 14 th minutes 
of the simulation. The amount of these increases 
is determined solely by the model validation algo- 
rithm. A final level of 8% is required to validate the 
model at the flight condition of 18 minutes. This 
level agrees with Figure 8 which shows a variation 
in damping between the nominal and true plants 
of approximately 8% at this flight condition. 

The modal validation portion of the flutter analy- 
sis was stopped due to low damping values of less 
than .005 after 19 minutes. The algorithm encoun- 
ters some conditioning problems when the system 
is extremely close to an instability. 


8. Simulation : On-Line Flutter Margins 

Flutter margins are computed at each test point 
based on the levels of uncertainty required to 
validate the model with flight data from that 
test point. Figure 6 demonstrates these calculat- 
ed robust and nominal flutter margins along with 
the true flutter margin. 



Figure 6. Flutter Margins expressed as Dynamic Pressure 
Flight Condition : Robust /i Prediction ( — ), Nominal p Pre- 
diction (. . .), TVue Flutter Margin ( ) 


The robust flutter margin is not computed with 
respect to any uncertainty for the first 3 minutes 
so the nominal and robust margins are identically 
poor at first. The nominal margin does not account 
for any uncertainty and remains constant through- 
out the flight. The robust margin varies with time 
due to changes in the uncertainty levels required 
to validate the flight data at each test point. 

The robust margin decreases to remain conserva- 
tive to the true flutter margin due to the worst- 
case nature of the computation with respect to the 
increasing uncertainty levels. Note it is only conser- 
vative to the flutter margin for Ptrue at that partic- 
ular test point since the uncertainty set at time t\ 
does not necessarily cover the dynamics at a later 
time fa due to the time-varying mass variation. 

The /i-analysis method of computing robust flutter 
margins is able to utilize flight data from a stable 
test point. The robust flutter margin immediately 
provides a realistic flutter margin without requiring 
the aircraft to approach an unstable flight condi- 
tion. 

The nominal flutter margin is shown on Figure 6 
to demonstrate the large error in this flutter mar- 
gin as compared to the final true flutter condition. 
This nominal margin corresponds to the pre-flight 
prediction computed via /i or p-k methods with no 
accounting for modeling uncertainty. 
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The concept of a flutterometer is introduced as a 
tool for indicating the proximity to a flutter con- 
dition. One implementation of this concept is to 
describe the smallest increase in ^ which results in 
flutter. Figure 7 demonstrates this flutterometer 
concept for the simulated flight test. 



Figure 7. Flutterometer showing Distance to Flutter in 
Dynamic Pressure : Robust /i Prediction ( — ), Nominal /i 
Prediction (...), True Distance to Flutter ( ) 


The benefits of the based flutterometer are easily 
seen when considering the traditional method of 
predicting flutter margins based on the evolution 
of modal damping as shown in Figure 8. 



Tirm (min) 


Figure 8 . Modal Damping for the Trailing Edge Flap Mode : 
Nominal Model P n0 m (...), True Aircraft ( ) 

Tracking modal damping is clearly inferior to the 
/i method for on-line flutter prediction. Modal 
damping provides little information until minute 
17 when the aircraft is approaching flutter. After 
that time, damping shows a decreasing trend to 
indicate proximity to an instability; however, the 
nonlinear behavior precludes calculating the exact 
proximity of that instability. 


0. Simulation : Computational Cost 

The efficiency of any flight flutter test program is 
directly measured by the ability to compute accu- 
rate flutter margins in a minimal amount of flight 
time without sacrificing safety of the pilot and 
aircraft. This simulation demonstrates computing 
robust flutter margins can be efficiently performed 
in an on-line manner for flight test programs. 

The robust flutter margin procedure is imple- 
mented in Mat lab using the /i- Analysis and 
Synthesis Toolbox [1]. The computation times 
presented are for a standard 200 MHz Pentium 
computer often available in a flight analysis facility. 

The computation time at each test point results 
from validating the uncertainty description and 
computing robust flutter margins. These computar 
tion times are shown in Figure 9 with circles rep- 
resenting the analysis time at a specific test point. 



Figure 9 . Computational CPU Time to Determine Robust 
Flutter Margin at each Test Point (o) 

The 3 rd minute test point marks the first robust 
flutter margin computation. This initial uncertain- 
ty set is conservative to the time- varying dynamics 
until the 8^ minute, as evidenced by the constant 
computational cost. Similarly, the uncertainty set 
at the 14* h minute is conservative until the 18 th 
minute. The computational cost is reduced after 
this time because model validation is no longer per- 
formed due to the onset of flutter instability. 

The 8,10 and 13 minute test points require longer 
time than other test points. Figure 5 shows the 
uncertainty at these three test points is increased 
to account for the time- varying dynamics of Pt ruc . 
The additional analysis time corresponds to the 
extra iterations required for invalidating the old 
uncertainty levels and computing new levels which 
are not invalidated by the flight data. 
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The flight path shown in Figure 3 combined with 
Figure 9 is indicative of the efficiency of the robust 
algorithm; however, there are several aspects to a 
real flight test program that are not easily modeled. 
Firstly, the simulation did not take into account 
data transfer time between the aircraft and the 
analysis engines [6]. Also, the simulation assumed 
a single data set generated at each test point could 
be used for model validation while typically sever- 
al sets are required to generate a rich flight data 
set [9]. These issues are demonstrative of inefficien- 
cies associated with flight data which are being 
addressed through improved flight test procedures 
and distributed analysis environments [11]. 

Given a good data set, the robust flutter margins 
are computed for the majority of test points in less 
than 1 minute with the largest analysis time being 
only 2.5 minutes. This clearly does not present an 
excessive computational burden to the flight pro- 
gram. In fact, the benefits of the procedure are 
more emphasized when considering the analysis 
time to generate this accurate predictor is compar 
rable to the computational time required for tra- 
ditional damping estimates which provide far less 
information about the true flutter boundary. 

10. Conclusion 

Worst-case flutter margins computed with a robust 
stability method can provide valuable information 
about the flutter dynamics of an aircraft. This 
paper introduces an on-line implementation of such 
a method for use in a flight flutter test program. 

A simulated flight test demonstrates several advan- 
tages to the robust flutter margin method as com- 
pared to traditional flutter analysis methods. The 
robust method uses flight data from a stable flight 
condition to compute a flutter boundary without 
requiring the aircraft to approach the instability. 
Also, this method computes margins which remain 
conservative to time-varying dynamics by updat- 
ing uncertainty levels at every test point so the 
new flight data does not invalidate the model. 

The robust flutter margins quantify the proximity 
to a flutter condition which can not be comput- 
ed by tracking damping trends. This fact leads to 
the first formulation of a robust flutterometer tool. 
This method could dramatically increase efficien- 
cy of a flight test program by accurately predict- 
ing flutter margins despite unmodeled time- varying 
dynamics of the aircraft while increasing safety due 
to the conservativeness of the predictions. 
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